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Identification of a high- velocity compact nebular filament 
2.2 arcsec south of the Galactic Centre 
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ABSTRACT 

The central parsec of the Milky Way is a very special region of our Galaxy; it contains 
the supermassive black hole associated with Sgr A* as well as a significant number 
of early-type stars and a complex structure of streamers of neutral and ionized gas, 
within two parsccs from the centre, representing a unique laboratory. We report the 
identification of a high velocity compact nebular filament 2.2 arcsec south of Sgr A*. 
The structure extends over 1 arcsec and presents a strong velocity gradient of 200 
km s~^ arcsec"^. The peak of maximum emission, seen in [Fe III] and He I lines, is 
located at da = +0.20 ± 0.06 arcsec and d5 = —2.20 ± 0.06 arcsec with respect to 
Sgr A*. This position is near the star IRS 33N. The velocity at the emission peak is 
Vr = —267 km s~^. The filament has a position angle of PA — 115° ± 10°, similar to 
that of the Bar and of the Eastern Arm at that position. The peak position is located 
0.7 arcsec north of the binary X-ray and radio transient CXOGX J174540. 0-290031, 
a low-mass X-ray binary with an orbital period of 7.9 hr. The [Fe III] line emission 
is strong in the filament and its vicinity. These lines are probably produced by shock 
heating but we cannot exclude some X-ray photoionization from the LMXB. Although 
we cannot rule out the idea of a compact nebular jet, we interpret this filament as a 
possible shock between the Northern and the Eastern Arm or between the Northern 
Arm and the mini-spiral "Bar" . 
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1 INTRODUCTION 

The central parsec of the Milky Way (MW) is a very special 
region of our Galaxy. It not only hosts the central supermas- 
sive bl ack hole, associated with Sgr A*, with w ell determined 
mass (|Ghez et al.ll2008l : [GiUessen et al.ll2009h . but also con- 
tains a significant number of early-type stars, including some 
in the Wolf-Rayet stage. 

The central cavity has a mini-spiral structure of 
ionized gas streamers, encircled by a circum-nuclear 
disk (CND) of molecular gas. Most of the ionized gas 
emission appears in Hydrogen and Helium recombination 
lines, very likely originating from ionization produced 
by the many young and hot OB stars; the equivalent 
to some ^ 250 09 stars would be required to a c count 
for the emission (ILacv. Achtermann fc Serabvnl 199ll: 
Shields fc Ferlandl [l994l : [Roberts. Yusef-Zadeh fc Goss 



199e 



20041 ) 



: lMorris fc Serabvnlll996l : |Paumard. Maillard fc Morris 
An anomalous structure, called mini-cavity, was iden- 
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tified in the mini-spiral, with a diameter of 2", about 3'!5 
to the southwest of Sgr A* (JYusef-Zadeh. Morris fc Ekera 
1 19891 . 1 19901 : |Zhao ct al.' Il99ll), as s ociated with strong [Fe 
HI] emission (,Eckart c t al.l 1 19921 : iLutz. Krabbe fc Genzell 
Il993h . For a recent review of the literature on th e Galactic 
centre, see lGenzel. Eisenhauer fc Gillessenl (|2010l ). 

X-ray binary t ransients are overab undant in the central 
parsec of the MW. lMuno et al.l (1200581 ) have found that the 
overabundance is at least a factor of 20, per mass, in the 
central parsec when compared to the abundance at a dis- 
tance between 1 and 23 parsecs. It is likely that this over- 
abundance is due to three-body interactions between binary 
star systems and either black holes or neutron stars that 
have concentrated in the central parsec through dynamical 
friction. This is also supposed to be the case with the for- 
mation of similar s ystems in the cores of globular clusters 
iMunoet al.ll2005al ). 

The known X-ray transient closest to Sgr A*, at only 
2"8 south of the centre of the MW, is CXOG X J174540.0- 
290031, discovered in 2004 (|Muno et al.ll2005al ). It was found 
to be an X-ray binary with an orbital period of 7.9 hr 
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IjMuno et al.ll2005al lb^■ From K band observations, an up- 
per limit for any counterpart was set as if < 16, suggesting 
that it is a low-mass X-ray binary ( LMXB). This X-ray tran- 
sient also revealed a radio outburst (jBower et al.ll2005l ) . This 
large radio outburst, as well as its X-ray emission, suggests 
that this source mo re likely contains a black hole than a neu- 
tron star primary (JFender fc Kuulkersll200ll : IPorquet et alj 
l2005h . 

In this paper we report the discovery of a high veloc- 
ity compact nebular filament extending ~ 1" with a strong 
velocity gradient that could possibly be associated with a 
colliding shock between the streamers of the mini-spiral. 



2 OBSERVATIONS AND DATA REDUCTION 

This paper analyzes archive data obtained with the Spectro- 
graph for INtegral Field Observations in the Near Infrared 
(SINFONI) on the VLT, available in a public access data 
ban4j. The observations were made in K band on April 22, 
2007. A fore-optics of 0"250 pixel"^ was used, providing a 
field of view (FOV) of 8" x 8". A total of four exposures 
of the Galactic nucleus were retrieved from the public data 
bank to be used in this work. 

Standard calibration images of linearity lamp, distor- 
tion fibre (used in the data reduction to compute spatial 
distortions and perform a spatial rectification), flat lamp, 
arc lamp and sky field were obtained from the public data 
bank. Calibration images of an AOV standard star were also 
retrieved from the data bank. The data reduction was made 
with Gasgano software and included the following steps: 
determination of the trim, sky subtraction, correction of 
bad pixels, flat-field correction, spatial rectification (includ- 
ing correction for spatial distortions) , wavelength calibration 
and data cube construction. At the end of the process, four 
data cubes were obtained, with spaxels (spatial pixels) of 
0"125 X 0"125. Afterwards, with a script developed by us 
and using data from the AOV standard star, we applied, in 
an IRAF environment, telluric absorption removal and flux 
calibration to all the data cubes. Details of the data treat- 
ment are given in Appendix A. 



3 DATA ANALYSIS 

After the data treatment, we applied PGA Tomography 
(see Appendix A) again. Tomogram 1, shown in green in 
Figure 1-a, displays the young stellar population. Tomo- 
gram 2, shown in blue in Figure 1-a, reveals stars with 
strong i nfrared excess , pres umably very young stars, like 
IRS 21 (JTanner et al.l 120021 ') and IRS 1. The red compo- 
nent in Figure 1-a corresponds to an image of stars obtained 
by subtracting the edges of the CO bands from the recon- 
structed data cube, after suppressing eingenvectors 1 and 
2 (see Steiner et al. 2009 for a description of "feature sup- 
pression"). The stellar component revealed by Figure 1-a 
basically shows the known young as well as the old stellar 
populations. IRS 16, the closest IRS source to the Galactic 
centre, is resolved in the individual stars, as is IRS 33 (SW 
of the FOV). 
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We performed spectral synthesis on the spectrum 
of each spaxel of th e datacube analyzed here , using 
the Starlight software (|Cid Fernandes et al.l [2005|), which 
matches the stellar spectrum of a given object with a model 
corresponding to a combination of template stellar spectra 
from a pre-established base. In this work, we used the base 
of stellar spectra MI LES (Medium resolution INT Library of 
Empirical Spectra) (jSanchez-Blazguez et al.ll2006l ). contain- 
ing the spectra of 150 stellar populations with ages between 
1.0 X 10® years and 1.8 x 10^'' years and with metallicities 
between 0.0001 and 0.05 (with Zq = 0.02). The spectral 
synthesis resulted in a synthetic stellar spectrum for each 
spajcel of the data cube. 

Before performing spectral synthesis on the data cube, 
however, we prepared the spectra in two steps: correcting the 
interstellar extincti on due to the Galaxy, using Ay — 1 9 and 
the reddening law of lCardelli. Clayton fc MathisI (|l989l ) ; and 
carrying out a spectral re-sampling with AA = lA per pixel. 
The first step used a script written in IDL language, while 
the second used the task "dispcor" , of the "noao" package, 
in IRAF environment. 

In order to analyze the emission line spectrum in more 
detail, we subtracted the synthetic stellar spectra obtained 
with the spectral synthesis from the observed ones. This 
procedure resulted in a data cube with emission lines only. 
We then constructed RGB images of the emission lines Br7 
and He I A20581j4 of this emission line data cube. We can 
see that the recombination lines Br7 and He I A20581j4 are 
detected as stellar winds in early-type stars and, also, as 
ionized interstellar lines basically delineating the Northern 
Arm. A number of other lines are also seen in the southwest 
of the FOV: A2145lA, A22178A, A2242lA, and A23479A are 
the strongest among them and are identified as [Fe III] lines 
(Lutz et al. 1993) . We then constructed RGB images of the 
emission lines Br7 and [Fe HI] A22180A, shown in Figure 
1-c and Figure 1-d, respectively. 

We identified a high velocity compact nebular filament, 
2"2 south of Sgr A*. The structure extends over ~ 1" and 
presents a velocity gradient of ~ 200 km s^^ arcsec"^ for He 
I, Br7 and [Fe HI] A22180A lines (Figure 1-e, and Figures [H 
and \^ . The peak of maximum emission seen in the [Fe III] 
and He I lines is located at da = +0"20 ± 0"06 and dS = 
-2"20 ± 0"06 with respect to Sgr A*. This position is near 
the star IRS 33N, at da = -0"06; dS = -2"19 from Sgr 
A*. As can be seen in Figure |3] the velocity at the peak, as 
estimated from the Br7, He I A2058lA and [Fe HI] A22178A 
is ^ —267 km s~^. The filament has a position angle of PA 
= 115° ± 10° (see Figures [2] to H]) . 

The observed line ratios at the emission peak of the 
filament are (He I A2058lA)/Br7 = 0.39 ± 0.06 and ([Felll] 
A22180A)/Br7 = 0.27 ± 0.04. 



4 DISCUSSION AND CONCLUSION 

The filament is positioned in one of the mini-spiral regions 
nearest to the Galactic centre, in the field where the East- 
ern and the Northern Arms of the mini-spiral cross. It is 
also at the eastern edge of the mini-cavity. This field shows 
significant blue-shifted He I emission line (Paumard et al . 
2004) and also in the H92a and H30a lines (Zhao ct all2009l . 
l2010l i. 
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Figure 1. a) The FOV of the SINFONI observation in pseudo colors. In green we show the young stars; in red, the old stars and in blue, 
the stars that have strong infrared excess, b) Contours of the stars in the FOV, with identifications, c) The FOV of the Br7 emission line 
in three colors, showing the radial velocity structure of the Northern Arm. The colors blue, green and red represent the velocity ranges 
-156 km s"^ ^ K ^ -73 km s'^, -46 km s"! s£ K sj -18 km s'^ and 10 km s^^ s£ V,- ^ 93 km s"\ respectively d) The same for 
[Fe III] A22I8OA. The colors blue, green and red represent the velocity ranges —108 km s^^ Sj Vr ^ km s^^, 27 km s^^ 55 Vr Sj 54 km 
s~^ and 81 km s~^ ^ V^ ^ 189 km s~^, respectively, e) The filament, seen in the line [Fe III] A22180A. The colors blue, green and red 
represent the velocity ranges -379 km s^^ < Vr ^ -298 km s-\ -271 km s'^ ^ K < -243 km s"! and -216 km s'^ ^ V^- ^ -135 
km s"-'^, respectively, f) Contours of the image shown in e). The + sign indicates the position of Sgr A*. 



The filament that we see may be described as a "jet- 
like" structure, being brighter at the centre and fading to- 
wards the borders, at least as it is seen in the He I and 
the [Fe III] lines (see Figure [3}. This is in the part of the 
mini-spiral closest to Sgr A*. The challenge is to explain the 
compactness of the structure as well as its ionization (He I) 
and excitation ([Fe HI]). 

The filament overlaps spatially with the star IRS33N 
(see Figure 1-b and Figure |4ll. Could this star be the origin 
of a jet? IRS 33N is a if = 11.1 star, presenting a spectral 



type of BO. 5-1 I and has a rad ial velocity of Vr = +68 ± 20 
km s~^ (|Paumard et al.l I2OO6I ) while the centre of the fil- 
ament is at Vr = —267 km s~^. Although close in terms 
of spatial projected distance, the nebular filament is kine- 
matically decoupled from the star IRS 33N. Therefore this 
association can be discarded. 
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Figure 2. Chanel maps for Br7. The velocity of each channel is given on the top. The position of Sgr A* is marked as a + sign. The 
contours of Figure 1-f are shown in green. We call the attention to the fact that at highest velocities, the filament is the strongest emission 
detected in the FOV. 
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Figure 3. Left: Pseudo-slit of the data cube along the gaseous filament, displaying spatial position versus velocity for the [Fo III] A22180A 
line. Centre: The same for the line He I A20581A. The star at position —1" and velocity ~ +250 km s~^ is IRS 33SE. Note the presence 
of a second nebular velocity with Vr ~ +80 km s"'^ at the same position as the high velocity filament. This is the velocity of the Eastern 
Arm at that position. Right: The same for Br7. 
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Figure 4. The field of the nebular filaments shown at positions 
with distinct velocities: 1 (—378 km s~^); 2 (—323 km s~^); 3 
(-267 km s-l); 4 (-212 km s"!) and 5 (-156 km s"!). The 
northern edges of the Northern and Eastern Arm streamers, as 
well as of the mini-cavity, are shown. Three stars of IRS 16 as well 
as of IRS 33 are also shown. The position of Sgr A* is marked as a 
-|- sign. The error circle of the LMXB (green) and the associated 
radio transients (red) are also shown. 



4.1 FireAvorks in the mini-spiral: a shock between 
streamers? 

The orbits of the mini-spiral have been described as contain- 
ing a one-armed spiral in a Keplerian disk (Lacy et al. 1991) 
that properly describes th e [Ne II] A12.8 fiv a. An alterna- 
tive description is given bv lZhao et al.l ((20091), in which the 
whole mini-spiral is essentially described by three streamers 
in elliptical orbits around Sgr A*, on the basis of H92a ob- 
servations. The Northern and Eastern Arm streamers have 
high eccentricities (e — 0.83 ± 0.10 and 0.82 ± 0.15, respec- 
tively) and cross near the periastron. The formal solution 
provides a projected crossing between the two orbits about 
4" to the so uth of the Galactic cent r e. Thi s solution has been 
disputed by llrons. Lacy fc Richten (|2012l ) who find that the 
one-armed spiral pattern may be a better description of the 
[Ne II] nebular emission. The observations of both H92q and 
the [Ne II] lines have good spectral resolution but poorer 
spatial resolution, of ~ l"25 to 2", than the data reported 
here. They could not have identified such a strong velocity 
gradient in a filament so compact. 

Perhaps the best explanation for the existence of this 
filament is that it is associated to a shock between streamers 
in the mini-spiral. In the one-arm pattern, the position of the 
filament is where the one-armed spiral ends; that could also 
be seen as the crossing point between the Northern Arm and 
the Bar. In the three ellipses model, the filament is near the 
position where the projected orbits of the Northern and the 
Eastern Arm streamers cross. It is perhaps not a coincidence 
that the position angle of the filament (PA — 115° ± 10°) 
is similar to that of the Eastern Arm and also of the Bar. 
Perhaps the most challenging aspect of this interpretation is 



how to explain the large velocity gradient of ~ 200 km s~^ 
arcsec^ ^. It probably has to do with the proximity of the 
supermassive black hole and its steep gravitational potential 
gradient. 

The blue arm of the filament points toward the mini- 
cavity. This structure was modeled as a bubble of ionized 
hot gas generated by a fast (^ 1000 km s^^) wind originated 
froin one or several sources, within a few arcseconds from the 
mini-cavity, and blown into the partially ionized neutral gas 
streamers orbiting the centre of the Galaxy (Eckart et al. 
1992; Wardle & Yusef-Zadeh 1992; Lutz et al. 1993; Meha, 
Coker & Yusef-Zadeh 1996). 

In fact, af ter computing the elliptical Keplerian orbits of 
the streamers, IZhao et al.l (|2009l ) suggested that the North- 
ern and Eastern Arms may collide in the "Bar" region, a 
few arcesconds south of Sgr A*. This is the region with the 
highest specific kinetic energy of both Northern and East- 
ern Arms ( s ee Fig ure 21 in Zhao et al. 2009). According to 
IZhao et al.l ((20091), the high specific kinetic energy density 
emphasizes the possibility that the mini-cavity was created 
by strong shocks resulting from the collisions of the two 
streamers. 



4.2 The X-ray transient and origin of the nebular 
lines 

The peak of the emission of the filament is located 0"7 north 
of the binary X-ray transient CXOGX J1745 40.0-290031, 
a LM XB with an orbital period of 7.9 hr ((Muno et al.l 
l2005al (a). This source is associated to o ne transient and on e 
persistent albeit variable, radio sources ((Bower et al.ll2005l '). 
The X-ray transient source CXOGX J174540.0-290031 was 
detected on July, 5-7, 2004 and is lo cated 2"8 south of Sgr 
A*, with an uncertainty of about 0"3 ((Muno et al.ll2005al (H). 
Observations made at 43 GHz a few months earlier, on 
March 28, 2004, detected two sources, named NE (north- 
east, the strongest one, reaching 93 mjy at 43 GHz) and 
SW (southwest) with a peak flux of 48 mJy. The coinci- 
dental variability of both radio a nd X-ray sources strongly 
suggests that they are associated. (Zhao et al.( ((20091 ) also re- 
port radio observations at 1.3 cm and found that the SW 
source is a persistent one, being detected before the transient 
phenomenon at the 1.1 ±0.1 mJy level. The NE component, 
not detected in previous observations, was detected after the 
transient. These authors interpret the NE component as a 
hot spot associated with the head of the jet ejected from the 
persistent SW radio-variable source. If so, then the stellar 
counterpart of the X-ray binary is probably near the star 
IRS 33SW, close to the X-ray error box (Figure [4|. The po- 
sition angle of the radio jet is PA ~ 60° , quite distinct from 
the position angle of the nebular fllament. As the filament is 
at ~ 0"7 from the X-ray source and has a distinct position 
angle, we may rule out any direct association between the 
X-ray binary and the origin of the nebular filament. 

The SINFONI observations (April 22, 2007) were made 
a few years after the radio detection (March 28, 2004) and 
also after the X-ray transients detections (July 5-7, 2004). 
The nebular filament is characterized by a spectrum that is 
quite distinct from other gas emitting regions in the FOV. 
Br7 and He I A20581A are the strongest lines. Other gas 
emitting regions in the FOV are probably photoionized by 
the young stellar population (such as in IRS 16) that pro- 
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duces UV photons capable of ionizing H and He (Genzel et 
al. 2010). Such UV photons, because of their low energy, do 
not heat the nebulae and some 250 09 equivalent sta rs are 
necessary to account for the emission (|Zhao et al.ll201CI ) . The 
vicinity of the X-ray source is the only region in the FOV 
that presents such strong forbidden lines. But shock waves 
are also efficient in heating gas and enhancing forbidden 
line emissions. The field of IRS 33 is the hottest in the mini- 
spira l , as m easured from the H30a and H92a (|Zhao et al.l 
l2009l . i2010l '). with a temperature of 13000 ± 3000 K and an 
electron density of ~ 10^ cm"''. The temperature of 13000 
K is enough to excite the [Fe III] emission. It is not enough 
to ionize the gas to th e Fe"*"^ phase. This requires a temper- 
ature of T > 18000 K (|Osterbrock fc FerlandbOOd '). But the 
beam size of the radio observations is ~ 2"; thus a compact 
structure with T > 18000 K could certainly exist. For col- 
lisional ionization of helium, the temperature must be even 
higher: T > 30000 K. We conclude that shock heating could 
both produce and excite the Fe^"*" ion, but not ionize He. 
Therefore, a combination of non-shock ionization mechanism 
plus shock heating seems to be required. 

Could the X-ray transient be, at least partially, respon- 
sible for photoionizing the filament? Hard X-rays are efH- 
cient in photo ionizing gas and heating it at the same time 
(|Halpern fc Stcincr 19831 '). enhancing the emission from low 
ionization nebular species. However, the ionization parame- 
ter U must be [/ > 10~^ in order to produce high ionization 
species and U ^ 10"'^ in order to produce low ionization 
species (jFerland fc Netzerlll983l ). In the present case, the 
transient X-ray luminosity was I/a, ~ 4 x lO'^* erg s~^. How- 
ever, at the time of the transient, diffuse X-ray emission was 
observed ~ 2" — 3" southwest of the source. This suggests 
that the real X-ray luminosity was L^ ^ 2 x 10'^^ erg s^^ 
(jMuno et al.ll2005bl ). The discrepant X-ray luminosity seen 
directly is probably due to the fact that the X-ray binary is 
seen nearly edge-on. 

Assuming that the x-ray luminosity is La; ~ 2 x 10"^^ erg 
s~^ we obtain an ionization parameter of [/ ~ 5 x 10~^ at the 
filament; this is too low to produce a significant ionization 
effect. We conclude that the X-ray binary has no ionization 
or heating effect on the filament unless the source radiates 
soft X-rays near the Eddington limit. This could not have 
been directly observed due to the strong interstellar absorp- 
tion and would produce U ~ 10^^ at the filament and could 
be responsible for the strength of the [Fe HI] emission and, 
perhaps, of He I. 

Further studies are required to elucidate the nature of 
this filament and the role of the LMXB in its neighborhood; 
it is important to observe the FOV with better spatial res- 
olution. 



4.3 Summary 

We have analyzed a data cube centered ~ 2" southeast of 
Srg A*, obtained with the SINFONI/VLT IFU spectrograph 
and studied the kinematics of the gaseous ionized emission 
in the FOV. Our main findings are: 

• We identified a compact nebular filament, 2"2 south of 
Sgr A*. The structure extends over ~ 1" and presents a 
velocity gradient of ~ 200 km s~^ arcsec"^. There is a peak 
of maximum emission seen in [Fe HI] and He I lines, located 



at a; = -|-0"20 ± 0"06 and y = -2'f20 ± 0"06 with respect 
to Sgr A*. This position is located at ~ 0"3 east of the star 
IRS 33N. This emission has the most negative velocity in 
the FOV we have analyzed. 

• The velocity at the emission peak, as estimated from the 
Br7, He I A2058lA and [Fe HI] A22178A is ~ -267 km s^V 
The nebular filament has a position angle of PA = 115° ±10°, 
similar to that of the Eastern Arm at that position and to 
the position angle of the mini-spiral Bar. Although close, 
the nebular filament is, kinematically decoupled from the 
star IRS 33N, which has a velocity of Vr ~ +68 ± 20 km 

• The peak position is also located 0''7 north of the binary 
X-ray transient CXOGX J174540.0-290031, a LMXB with 
an orbital period of 7.9 hr. This source is associated with 
one transient plus one persistent and variable radio sources. 

• The [Fe HI] line emission is strong in the filament and its 
vicinity. These lines are probably excited by shock heating, 
lonization/heating by soft X-ray photons from the nearby 
LMXB is possible, although unlikely, because of the high 
luminosity required. 

• Although we cannot rule out the possibility that we 
are dealing with a jet-like structure, we propose that the 
filament may represent a shock from the collision between 
the Northern and the Eastern Arm streamers around the 
centre, or perhaps, between the Northern Arm and the Bar. 
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APPENDIX A: DATA TREATMENT 

After the data reduction, a correction of the differential 
atmospheric refraction was applied to all data cubes, using 
an algorithm developed by our group. It is important to 
mention, however, that the relative shifts observed in the 
images of all data cubes are not compatible with what is 
expected for the differential atmospheric refraction effect. 
Therefore, we believe that these relative shifts are, at least 
partially, a consequence of some instrumental effect and not 
entirely due to the differential atmospheric refraction. In 
order to combine into one the four data cubes obtained after 
the correction of the atmospheric differential refraction, we 
divided the data cubes in two groups with three data cubes 
each. Naming the data cubes 1, 2, 3 and 4 (according to 
the order of the observations), the division in two groups 
was the following: 



of the data cube, we obtained spaxels of 0"0625 x 0''0625. 
This procedure does not change the spatial resolution of the 
observation, but sharpens the contours of the structures in 
the images. The disadvantage of such spatial re-sampling is 
that it introduces high spatial frequency components in the 
images (usually in the form of small vertical and horizon- 
tal stripes). The high spatial frequency components, how- 
ever, are not a serious problem because, after the spatial 
re-sampling, a Butterworth spatial filtering was applied to 
all the images of the data cube, to remove high spatial fre- 
quency noises, including the high frequency components in- 
troduced by the spatial re-sampling. The steps of this pro- 
cess were as follows: calculation of the Fourier transforms 
{F{u, v)) of all the images of the data cube; multiplication of 
the Fourier transforms F{u, v) by the image corresponding 
to the Butterworth filter {H{u, v)); calculation of the inverse 
Fourier transforms of all the products F{u, v) x H{u, v); ex- 
traction of the real part of the calculated inverse Fourier 
transforms. The Butterworth filter used in this case corre- 
sponds to the product of two identical circular filters, with 
orders n = 2 and cutoff frequencies along the horizontal and 
vertical axis equal to 0.27 Ny (Nyquist frequency). 

We detected the presence of a probable instrumental 
"fingerprint" in the data cube, after the Butterworth spatial 
filtering, which appeared in the data cube as a large horizon- 
tal strip in the images and also had a low frequency spectral 
signature. In ord er to remove it, we a pplied the PGA Tomog- 
raphy technique (jSteiner et al.ll2009r ) to the data cube. PGA 
transforms data expressed originally in correlated coordi- 
nates into a new system of uncorrelated coordinates (eigen- 
vectors) ordered by principal components of decreasing vari- 
ance. PGA Tomography consists in applying PGA to data 
cubes. In this case, the variables correspond to the spec- 
tral pixels of the data cube and the observables correspond 
to the spaxels of the data cube. Since the eigenvectors are 
obtained as a function of the wavelength, their shape is sim- 
ilar to spectra and, therefore, we call them eigenspectra. 
On the other hand, since the observables are spaxels, their 
projections on the eigenvectors are images, which we call 
tomograms. The simultaneous analysis of eigenspectra and 
tomograms allows one to obtain information that, otherwise, 
would possibly be harder to detect. Using PGA Tomogra- 
phy, we were able to identify and remove the instrumental 
fingerprint of the data cube analyzed here (Steiner et al., in 
preparation) . 

Finally, a R ichardson-Lucy deconvolution (|Richardsonl 
1 19721 : ILucvII1974| ) was applied to all the images of the data 
cube, using a synthetic Gaussian PSF. The use of a Gaussian 
PSF was possible in this case because we verified that the 
original PSF of this observation had an approximate Gaus- 
sian shape. The final data cube obtained after the full data 
treatment has a point-spread function (PSF) with FWHM 
~ 0"56 and a spectral resolution oi R ^ 5875. 



Group 1: 1, 2 and 3 
Group 2: 2, 3 and 4 

We then calculated a median of the data cubes of 
each group, which resulted in two data cubes at the end of 
the process. Finally, we calculated the average of the two 
data cubes and obtained the final, combined data cube. 

By performing a spatial re-sampling to all the images 
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